HOMOGENIZATION OF SPECTRAL PROBLEMS 
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Abstract. Homogenization of a spectral problem in a bounded domain with 
a high contrast in both stiffness and density is considered. For a special critical 
scaling, two-scale asymptotic expansions for eigenvalues and eigenfunctions are 
constructed. Two-scale limit equations are derived and relate to certain non- 
standard self-adjoint operators. In particular they explicitly display the first 
two terms in the asymptotic expansion for the eigenvalues, with a surprising 
bound for the error of order e 5 / 4 proved. 
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1. Introduction 

Homogenization for problems with physical properties which are not only highly 
oscillatory but also highly heterogeneous has long been documented to display un- 
usual effects, for example the memory effects observed by E. Ya. Khruslov [9l ll3[[T4] . 
Of particular interest in this context are the double-porosity models where the pa- 
rameter of high-contrast S is critically scaled again the periodicity size e, S ~ e 2 , 
e.g. [2j H]. Those have been treated both by a high-contrast version of the classi- 
cal method of asymptotic expansions, e.g. [TBI El H2] and using the techniques 
of two-scale convergence, e.g. [TH1 HTJ1 IS]. In particular, for spectral problems in 
bounded [19] and unbounded [20] periodic domains V.V. Zhikov studied the spec- 
tral convergence, introduced two-scale limit operator, developed the techniques of 
two-scale resolvent convergence and two-scale compactness. In [12] the spectral 
convergence of eigenvalues in the gaps of Floquet-Bloch spectrum due to defects in 
double-porosity type media were studied, and [5] supplemented this by the analysis 
of eigenfunction convergence based on an analysis of a uniform exponential decay. 

In this work we study spectral problems of double-porosity type in a bounded 
domain f2 where the high contrast might occur not only in the "stiffness" coefficient 
but also in the "density" , and argue that this leads to some interesting new effects. 
Namely, referring to the next section for precise technical formulations, for the 
spectral problem 

(1) - div (a £ (x) Vu £ ) = X £ p £ (x)u e , 

with Dirichlct boundary conditions on the exterior boundary, most generally, both 
a s and p e are e-periodic, a £ = p e = 1 in the connected matrix and a £ ~ e Q , p e ~ 
e 13 in the disconnected inclusions. (Outside homogenization, the above resembles 
problems of vibrations with high contrasts in both density and stiffness, e.g. [3]) 
The double-porosity corresponds to a = 2 and (3 = 0. For /? ^ 0, it is not hard to 
see that it is a = f3 + 2 when the spectral problems at the macro and micro-scales 
are coupled in a non-trivial way. To explore this, we choose = — 1 and a = 1 
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Figure 1. The geometry and the periodicity cell 



and show that this leads to some unusually coupled two-scale limit behaviors of the 
eigenfunctions and the eigenvalues. 

Namely, although the limit behavior of the eigenfunctions is still somewhat sim- 
ilar to that of double porosity, i.e. the two-scale limit is a function of only slow 
variable x in the matrix and a function of both x and the fast variable y in the 
inclusions, the limit equations themselves are quite different. We show that there 
exist asymptotic scries of eigenvalues A £ ~ Xq + e\i with Ao being any eigenvalue of 
a non-standard self-adjoint "microscopic" inclusion problem, Theorem 13. 1[ whose 
eigenfunctions are directly related to the two-scale limit wo(x,y) in the matrix. In 
fact, Ao is cither a solution of (3(\q) = |Qi|Ao, where /3(A) is a function introduced 
by Zhikov [19] , or is an eigenvalue of the Dirichlet Laplacian in the inclusion Qo with 
a zero mean eigenf unction. In the matrix, u e ~ v (x), where v is an eigenfunction 
of the homogenized operator in fi, whose eigenvalue v determines the second term 
Ai in the asymptotics of A e , see (|57p . This is first derived via formal asymptotic 
expansions, but then we prove a non-standard error bound: 

|A £ -A -eAi| < Ce 5/4 , 

see Theorems 14.61 & 14.71 The proof employs a combination of a high contrast 
boundary layer analysis with maximum principle and estimates in Hilbert spaces 
with £-dependent weights. We finally briefly discuss further refinement of the results 
via the technique of two-scale convergence. Namely, some version of the compact- 
ness result holds, cf. [19], indicating at the presence of gaps in the spectrum for 
small enough e, see Theorem 15. II 

The paper is organized as follows. The next section formulates the problem and 
introduces necessary notation, Section 3 executes formal asymptotic expansion and 
derives associated homogenized equations. Section 4 proves the error bounds and 
Section 5 discusses the two-scale convergence approach. Some technical details are 
assembled in the appendices. 

2. Problem statement and notations 

We consider a model of eigenvibrations for a body occupying a bounded domain 
in R™ (n = 2, 3, . . . ) containing a periodic array of small inclusions, see FigureQ] 

The size of inclusions is controlled by a small positive parameter e, e — > 0. First we 

introduce necessary notation. 

Let Q — [0,1]™ be a reference periodicity cell in W l . Let Qo be a periodic 

set of "inclusions", i.e. Qo + m = Q Q , Vra S Z™, and Qo — Qo f~l Q is a reference 

inclusion lying inside Q with C 2 -smooth boundary T, see FigureHJ Let Q\ — Q\Qo, 
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Qi = M. n \Q , T = dQo = dQi. Introducing y = x/e we refer to y as to a fast 
variable, as opposes to the slow variable x. In the x-variable the periodicity cell is 
eQ = [0, e) n . If y G Qj then i = eye eQj, j = 0, 1. We denote Slg := SI n eQ , 
Qf := Sin eQi = Q\nf, T e := ef n SI, see Figure HJ The trace on T E of function 
/ : Clj — > M™ is denoted by Let rtj, be the outer unit normal to Qq on its 
boundary T and let n x denote the similar normal on T E . 
Let stiffness a £ and density p £ be as follows 

, . / i, x e sif . / l, x e sif 

v ; \ e, x G Qg F£W e , x G S1q 

with a small positive e. 

We study the asymptotic behaviour of self-adjoint spectral problem 

(2) / a e (x)Vu £ V(j)dx - \ £ p £ (x) u £ <\> dx = 0, V</> G i?o(Sl) 

as e — > 0. If r and <9S1 are smooth enough then variational problem @ can be 
equivalently represented in a classical formulation 

(3) — div (a e (x) Vu £ ) = X e p £ (x)u £ , x G SI, 

(4) w £ |ao = 0, 

implying that at the interfaces the transmission conditions are satisfied 

du £ 



(5) 



l 



3. Formal asymptotic expansions 

We seek formal asymptotic expansions for the eigenvalues A 6 and eigenfunctions 
u £ in the form 

(6) A e ~ A +£A 1 +e 2 A 2 + ..., 



(7) u £ {x) 



v \ x, + svi (x, + s 2 v 2 yx, + . . . , a; G Of, 
w (x, + eioi ^x, + e 2 w 2 ^ar, + . . . , x G Slg. 



Here all the functions Vj(x,y), Wj(x,y), j > 0, are required to be periodic in the 
"fast" variable y; vq and wq are not simultaneously identically zero 

(8) vl + wl^0. 

In a standard way, the ansatz ©, (J7J) is then formally substituted into 
In particular, from ([3]), for (x,y) 6 fl X Qi, we obtain 

(9) -A y v = 0, 

(10) A„r, = 2 t ~''" 



9^ 9yj ' 



<9 2 i>i 

(11) -AyV2 = 2 Q x .Q y . + ^Z^O + ^0«0, 

(with A a and A x denoting the Laplace operators in y and x, respectively, and 
summation henceforth implied with respect to repeated indices), and for (x,y) G 
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f2 x Qq we have 

(12) - A y w 

(13) -A yWl = 2 



(14) 



-AyW 2 



X w , 

dx.jd yj 
d 2 w\ 



dxj dyj 

Further, the first of conditions ([5]) transforms to 



+ Aiw + Ao^i, 

+ A x wq + A 2 wo + Aiwi + \qW2- 



(15) 



Vj(x,y) =Wj(x,y) 
yer 



yer 



, x G n, j = 0,l 



Similarly, the other transmission condition ([5]) yield 
(16) 

(17) 

(18) 



dv 
dn y 


y er 


= o, 








dv\ 




dv 




dw 




dn y 


y er 


dn x 


H 

y er 


dn y 


yer 1 






dvi 




dwi 




dn y 


y er 


dn x 


H 

y er 


dn y 


y er 



dw 



y er 



The above has employed the identity 



(19) 



where 



du ( x\ i Uu t)U . . 



_i du 



du 

dn x 



x 

V = -, 

e 



i y ■ V a , g^- := n y ■ V X) with V a and V K standing for gradients in y 



and x, respectively. 
Finally, (j4]) suggests 

(20) 



xedn 



w 



xedn 



0. 



(The boundary layer problem does not generally permit satisfying Q by u, and uij 
for j > 1, as also clarified later.) 

Combining (|9|) and (|16p . together with the periodicity conditions in y, implies 
that vq is a constant with respect to y, i.e. 

Vo(x,y) = v°(x). 

Then, (|10|) and ()17j) form the following boundary value problem for v\ 
(21) -A y v 1 (x,y)=Q in SlxQj, 
The latter is solvable if and only if 



dv\ 
dn„ 



dv° 


H 


dw Q 




y er dn x 


y er 


dn y 


y er 



(22) 



r dn y 



dy = Q. 



Considering next (|12p and (TT5)) gives 
(23) - A y w = A w in fix Q , w (x,y) 

Since 



7r 9n y V J Ql 



A v wq dy = -A / w 

Qo 
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condition ([22]) is equivalent to 

(24) A ( Wo ) = 0, 

where 



(u) := / u(y)dy. 

Jqo 

We notice that ([2"3 |) -(f2"4 |) together with ([8|) constitutes restrictions on possible 
values of Ao- Those are described by Theorem 13.11 below. Before, let us consider an 
auxiliary Dirichlct problem 



(25) -A y( ^ = \ D <j) in Q , 4> 



0. 



Let {XjYjLx be eigenvalues for (|25[) . labelled in the ascending order counting for the 
multiplicities, and let {4>j}j^i be the corresponding eigenfunctions, orthonormal in 
L 2 (Qo), i.e. 



/ <P 3 (f>kdy — Sjk, 

J On 



where Sjk is Kronecker's delta. Denote by an the spectrum of (f2"S"|) : ao = \SjL\ 
We additionally introduce the following auxiliary problem: 



(26) - A. y r/ = A J? in Q , 7](y) 



= 1. 



Notice that ([26)1 is solvable if and only if Ao ao or Ao = with all the associated 
eigenfunctions <f>j having zero mean, (cj)j) — C0- In the former case 77 is determined 
uniquely and (|23[) implies wo(x,y) = v°(x)rj(y). In the latter case ij is determined 
up to an arbitrary eigenfunction (j>j associated with Aj 3 , however (77) is determined 
uniquely. 

By direct inspection, (|2"3"1) . (|2^|) has a non-trivial solution (i! ,!^), i.e. with ((SJ) 
holding, if and only if Ao is an eigenvalue of following problem: 

(27) -A. y C-AoC in Q , C(v) = constant, A (C) = 0. 

Theorem 3.1. The problem (1271) is equivalent to an eigenvalue problem for a self- 
adjoint operator in ^(Qo) with a compact resolvent. Therefore the spectrum of 
(|27l) is a countable set of real non-negative eigenvalues (of finite multiplicity) with 
the only accumulation point at +00, with the eigenfunctions complete in L2(Qo) 
and those corresponding to different Ao mutually orthogonal. 

The spectrum consists of all the eigenvalues X D of problem (1251) with a zero mean 
eigenfunction and all the solutions of the equation 

(28) B(Ao) := A (r?) = A f \Q \ + A f] Jf^ I = 



j'=i 



Xj -A 



(which are hence all real non-negative). In (|28[) the summation is with respect to 
only those Aj 3 for which there exists an eigenfunction with a non-zero mean. 

The associated eigenfunctions £ are either proportional to rj as in 
eigenfunctions of (|25p with zero mean. 



^We remark that the case of eigenvalues with zero mean is known to be not a "generic" case, 
i.e. unstable via a small perturbation of the shape of Qo, see e.g. discussion in 1101 and further 
references therein. 
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Proof. We claim that ([27]) corresponds to a self-adjoint operator associated with 
the (symmetric, closed, densely defined, bounded from below) Dirichlct form 

(29) a(C,h) := / VC • Vh dy 

JQo 

with domain 

(30) D(a) := {he ^(Qo) : h = constant}. 

To see this, in the weak formulation of the eigenvalue problem associated with 

(31) [ VC • V/i dy = A [ (h dy, Wi e D{a), 

JQo JQo 

we first set h to be an arbitrary function from C^°(Qo) which implies — A a £ = AoC m 
Qo, and then set h = 1 yielding Ao(C) = 0. Further, since the resolvent is obviously 
compact, each eigenvalue has a finite multiplicity, the set of all eigenfunctions £ is 
complete in L2(Qq) and those corresponding to different Ao are mutually orthogonal. 

Obviously, the spectrum of ([2"T]) includes those and only those eigenvalues of 
(1251) which have an cigenfunction cj>j with zero mean. In this case corresponding 
eigenfunctions of ([2"T]) are given by Q = C(j)j, 0. If Aj 5 does not have a zero- 
mean eigenfunction, then the solvability of (|2"T|) requires ( =0 implying £ = 0. 

yer 

Considering other possibilities, fix Ao outside ao and let rj be the unique solution 
of (|26]) . Then Ao is an eigenvalue of (|27|) if and only if 

(32) Ao(r/} - 0, 

with corresponding eigenfunction given by C,(y) = Crj(y), C ^ 0. 

Via the spectral decomposition, the solution to (|26[) is found to be, cf. [19] : 

_) 

3=1 "1 

Substituting ([33]) further into ([32]) yields ([28]). □ 
The formula ([2"5]) can be transformed to read 

(34) fl(Ao) =/3(A )- |Qi|A = 0, 
where function /3(A) has been introduced by Zhikov [19j : 

(35) /3(A) = A + A 2 ^ ( ^' )2 



(33) r 7 ( 2/ ) = l + A V-^-^(y). 

A ■ - A 



see Figure [2] This implies that Ao is either a solution to the nonlinear equation 
(36) /3(A) = |Qi|A, 

as visualized on Figure^ or is an eigenvalue of (|25[) with a zero mean eigenfunction. 

Remark 1. If Qo is a ball of radius < a < 1/2, i.e Qo = B a = {y : \y\ < a} + yo, 
then we have an explicit representation for /3(A). Indeed, for Ao o\d the solution 
of (|26|) is radially symmetric and (placing the origin in the ball's centre) reads 

rj(y) = |y|^ J^(A /2 |y|) (\a\^ J^X^a))^ , 
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where J n -2 (\y\) is Bessel function. Further, we have 



B(X ) = \ (v) = - 




= --\T\(l- n/2 + aXy 2 J',_ 1 _ 1 {\y 2 a)/J^_2{\y 2 a)) . 
Using (|35|) . (|33[) we obtain 

(3(\) = A(l - \B a \) - -|r| fl - n/2 + aX 1 ' 2 J'^ 2 {X 1 ' 2 a)/J^ 1 {X l ' 2 a)\ . 
/rt particular, for n = 3 we have, 

B(X ) = Xo(ri) = Ana (l — aX^ 2 cotan (A^ 2 a)^ , 

/3(A) = A(l - 47ra 3 /3) + Ana (l - aX 1/2 cotan (A 1/2 a)) . 

We next explore in detail the further steps in the method of asymptotic expan- 
sions, to determine v°, etc. Let us consider a i^-dimensional eigenspace (K > 1) for 
a given eigenvalue Ao of (|27|) . and let Ci, ■ ■ ■ , Ck" be associated linearly independent 
eigenf unctions. Then, (f2"3")) and (|2"1)) imply 

A" 

(37) woO,y) = 2j Cfe ^^ fe ^)- 

k=l 

Following Theorem 13. II we distinguish two cases: 

(a) A o\d. In this case (|26|) and (j23j) suggest 

(38) w (x,y) = v°(x)r)(y), 

and © implies «° ^ 0. 

(b) Ao G (t_d. The latter means Ao = Aj 3 for some j. This includes two further 
possibilities: 

(i) The eigenspace of (|25|) has an eigenfunction </>* with a non-zero mean. 
Since the solvability conditions for (|23|) include 

(39) «° = 0, 

necessarily v° = 0. Moreover, with Kr> denoting the multiplicity of Aj 3 
as of the eigenvalue of the Dirichlet problem (|25| . necessarily Kq > 2: 
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if Kd = 1 then wo = C(x)4>* and thus (j24|) implies C(x) = and wo 
contradicting to |[8]). Hence wo is given by (j37|) with K = Ku — 1, 
with fc = 1, if being linearly independent eigenfunctions of (|25|) 
with zero mean (such X eigenfunctions exist), 
(ii) All of the eigenfunctions corresponding \® have a zero mean. In this 
case wo is again given by (|37|) . with K = Ku if (??) i.e. B(\ ) ^ 
and if = if_D + 1 if -B(Ao) = with (k d +i = V where r)(y) is any 
solution of (|26|) . 



3.1. Case (a): Ao ^ ou- In this case Ao are solutions of There is a countable 
set of Ao = fij, j = 1, 2, . . . as Figure [3 illustrates. Note that this includes Ao = 0. 
Function /3 blows up at the points Aj 5 , which are eigenvalues of (j25|) having an 
eigenfunction with a non-zero mean, monotonically increasing between such points. 
It also directly follows from |35|) that /3(A) > |Qi|A for A G (0, Af), implying 
A? < ^2 < Af- Let Ao satisfying (|36|) be fixed. 

We consider problem (|2~lj) taking into account (|38|) . i.e. 



(40) -A yVl (x,y)=0 in fixft, 



yer 



dv 



yer 



0/ \ dr > 



yer 



where rj(y) solves (f2l))) and is given by ([55]) . Hence v\ is a solution to a problem 
depending linearly on v° and V^v , implying 



(41) 



dv° 

vi{x,y) = v°{x)N{y) + —Nj(y) + v$(x), 



with an arbitrary function Vi(x). The choice of v\ does not affect the subsequent 



constructions, so we set for simplicity v\ 
solutions to the problems 



0. In HTJ functions N* and Af are 



(42) 

and 
(43) 

Solvability of 



A y N j (y)=0 in Q u 

A y A%) = in Q u 
requires 



dN 
dn v 



-nj(y), 



drj 

yer dn y 



— dy 

r dn y 



0, 



which is equivalent to (f32|) and is hence already assured. Since the solutions of (J42 
and (I43p are unique up to an arbitrary constant, we fix those by choosing 



N^dy = / Af(y)dy = 0. 
Qi JQi 

We next consider the problem for wi, which from (|13p and (|15p combined with 
38l) reads 



(44) 
(45) 



-AyWi — XqWi = \\v rj + 2 



dv° dr/ 
dxj dyj 



in x Qo, 



= v°M(y) 



,er + ^ (y) 



yer 
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Since the problem depends linearly on v°, Xiv° and the solution admits repre- 

dv° 



sentation 
(46) 



wi{x,y) 



dx 



■{x)MM + v°(x)V(y) + \iv° (x)lZ(y) , 



where functions M.j, V and 1Z are solutions to the problems 



(47) 

(48) 
and 
(49) 



- A y M 3 - X Mj = 2p-{y) in Q , 
oyj 



A v V-X V = in Qq, V 



A y K - XqTZ = i](y) in Q , 



n 



o. 



Since by the assumption Ao ctd, all the problems (|47|) - (|49j) are uniquely solvable. 

The problem for V2 is in turn given by and ()18|) , whose solvability condition 
hence reads 

d 2 vi \ f ( dv\ dw\ dwo 



(50) 



A x v + X v + 2 



dxj dy-j 



dy = 



dn,. 



di 



dn. 



dy, 



with functions v\ , W\ and wq given by (|4ip , (|46|) and (|38| respectively. 

Appendix A provides a detailed calculation showing that the above yields the 
following equations for v°: 

(51) - divA hom V^° = i/(Ai)«° in ft, 

(52) 



on 



0. 



J i,k=i 



Here A hom = ^ h P m V 
rated domains, see e.g. [H] 

(53) A)r = \Qi\5jk 



is the classical homogenized matrix for periodic perfo- 



dNk 



dy; 



(54) 

where 
(55) 



i/(Ai) =CAi + A (|Qi|+ f Vdy), 



C I'- 



ll dy > 0. 



Qo 



Note that the problem (f5Tj) ~(f52" |) involves v — v{X\) as a spectral parameter. 
The spectrum of (f5"T]) -(f52" |) consists of a countable set of eigenvalues 

(56) < v\ < V2 < ■ ■ ■ < v n < > +oo. 

Corresponding eigenfunctions v n form an orthonormal basis in L2(f2), 

?;° ?;° dr — S 



Fixing an eigenvalue v of (|5ip , (|52p with corresponding eigenfunction i>° of unit 
norm in li%(££), according to (|54[) we find 



(57) 



A, = C" 



(v-Xo(\Qi\+J Vdy) 
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The following diagram summarizes the algorithm for constructing the first terms 
of the asymptotic expansions (for the case Ao $ <jd) 



J53J 

( f26t 



.4 



horn 



ED 
ED 



v, v 
Af 



o 



J57b 

Bit 

i47l-t49l 
l38l 



Ai 

M,v,n 

Wq 



RSI ITU, TO 

— > Wl — * l>2- 



We can additionally construct u>2 from (|14p and (115|) , whose unique solution exists 
for any choice of A2. For purposes of the justification of the first two terms in the 
asymptotics (the next section) it is sufficient to set A2 = and fix the corresponding 
solution W2- 

This completes constructing a formal asymptotic approximation, which we now 
summarize. We introduce an approximate eigenvalue 

(58) A £ = A +£Ai, 

and corresponding approximate eigenfunction 

v°(x) + evi (x, - ) + e 2 v 2 (x, — ) , 
f x \ ( x \ 2 £ ( 

Wq \ X, —j + £W\ \ X, — j + £ W2 \ X, 



(59) W s (x) 



x e fif , 
x e Qq. 



The essence of the above formal asymptotic construction is that the action of 
differential operator A £ on W £ defined by 



(60) 



A E W e := div (a £ VW E ) + A ePe W £ 



produces a small right-hand side in both fif and f2g, and on the interface T E in the 
following sense. 

Lemma 3.2. (i) max | div (a £ VW e ) + A e p e W e \ < Ce. 
(ii) max I div (a E VW E ) + A s p e W e \ < Ce 2 . 



(Hi) max 



dW £ 




<3W e 




on 





on 


1 



< Ce 2 



Proof, (i) Since the function W e is two-scale by the construction, in ilf 

div(a e vwy + A s/ 9 e w s = 

9 2 



e~ 2 A y +e~ 1 2 



+ A x + Ao +eAi^ (v°(a;) +e«i(x,j/) + e 2 v 2 (x,y))\ l 



(61) 



= fe- 1 A w «i(o;,y)+e f A y v 2 + + + X v°) + 

+ e 1 I 2 d " 2 — 1- A^ui + Ain° + A vi ) + e 2 (A^u 2 + Ai?;i + A0U2) + e 3 Aiu 2 1 

V 0^0% / J „= 



Since t>i is a solution to (|40| . the coefficient of e -1 vanishes. The same is with 
the coefficient of e° since V2 satisfies (|TTj) . Functions v°, V\ and V2 are solutions of 
elliptic problems with smooth enough coefficients to guarantee belonging solutions 
to C 2 . Thus, maxima for coefficients of e 1 , e 2 and e 3 in (lUTl) exist. 
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(it) Similarly, in Qq 

div (a e VWe) + AePsWe = 

= ^e _1 A y + 2 Qy + £ & x + e _1 Ao + Ai^ (wo(x, y) + ewi(x,y) + e 2 w 2 (x,y))\ y= ^ = 

= |e _1 (A y wo + X w ) + £° ( AyWi + 2 — % + XqWi + Xiwo } + 
[ \ dxjdyj J 

i ( a n d 2 wi . , \ 

+ e A y w 2 + 2- — h A x w Q + X w 2 + AiWi + 

V oxjdyj J 

+ e 2 (2 f y i 2 + A xWl +X lW2 ) + e 3 A x w 2 \ 

Since wo(x,y) = v°(x)r](y) is chosen according to (|2"6")) . the coefficient of e -1 van- 
ishes. The coefficient of e° vanishes due to (j4"4")l . Further, u>2 satisfies (|14[) with 
A 2 = and thus the coefficient of e 1 is zero as well. Since w\ and w% are solutions 
of elliptic problems with smooth enough coefficients, the maxima of the coefficients 
of e 2 and e 3 exist. 

(Hi) Using (HHJ), we obtain 



a, 



dn 



dn 



t^- +£ 'Q~~ S j ( w o( x 'V) +ewi(x,y) + e 2 w 2 (x,y))\ xere 



(s + -^~ \ (v° (x) + ev 1 (x,y) + e 2 v 2 (x,y)) 



\ On,, ihi,, <>n, J ; : -f s +£ {driy + dn x dn y dn x ) l*ff' + 



dvi 


dv°\ 




dn y 


dn x ) 


ver 


dwi 


dv 2 \ 




dn x 


dn x ) 


ier 
yer 



o / dvjQ dv\ dv° \ I j / 9wi 9wo dv 2 dv 

The coefficients of e° and e 1 vanish because of (|2"Tj) and (fT8|) respectively. The rest 
of the coefficients are smooth enough to guarantee that their maxima for x G T e 
and y G T exist. □ 

3.2. Case (b): A = Aj 3 . For simplicity, we consider here only the case of eigen- 
values of multiplicity K = 1 with zero mean eigenfunction (</> = <fij), assuming 
additionally Ao is not a solution of (j3"l| . All other degenerate cases, see page [JJ 
could be considered similarly. 

In this case we can introduce a refined approximation for the eigenfunction 

fevi (x, —] + e 2 v 2 (x, — J , x £ fif , 

ui ^x, —J + ew\ [x, -J+ew2 [x, — J , a; € S2§. 

where 

(64) w (x,y) = c(x)(j>(y). 

Lemma 3.3. Let c G C 3 (f2), i/ien £/iere ezisi smooth functions Vi,Wi, 1)%, W2 and a 
constant Xi such that A £ = X® + eXi and W* defined by (|63p satisfy 

(i) W*(x)&C(Q), 

(ii) max | div (a e VW e *) + A £ p £ W £ *| < Ce, 

Slf 

(in) max | div (a e VW*) + A e p E W*\ < Ce 2 , 
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fiv) max 



Proof. See Appendix B. 



dW* 




aw* 




£ dn 





an 


1 



<Ce 2 



□ 



4. Justification of asymptotics 



4.1. Operator formulation. We use a standard notation for Lebesgue and Sobolev 
spaces: Lp(Q,) is a p- weighted L 2 -space of square- integrable functions in CI. Notation 
(•, -)h is used for a scalar product in a Hubert space H. 

Let C £ — L 2 e (fi) and Tt £ be ifg (0) Sobolev space with a scalar product 



(u,v) 



/ a E (x)\7u ■ Vu dx + / p e (x)uvdx. 
Jn Jn 



Following a standard procedure, see e.g. [llj . we introduce a bounded operator 
B £ : C £ -»• £ e such that 

(65) (B e f,v) w = (f,v) c >, VveH £ . 

In other words £> e / = u e , where m £ is the solution of the problem 

(66) - div (a £ Vit e ) + p e u £ = p e f, x G CI, 

(67) u e \ m = 0, 

du E 

= U F , — - = £ 

1 



(68) 



dn. 



o' ] 

Note that operator B e is positive, self-adjoint and compact for any fixed e > (since 
its image is in Ti e ). Eigenvalue problem $2$ is equivalent to 

(69) B e u E = (A e + 1)"V in C e . 

Hence the spectrum of the problem consists of a countable set of eigenvalues 

< \\ < \ e 2 < ■ ■ ■ < \ E k < > +oo, 

with the only accumulation point at +oo. Moreover, the set of corresponding eigen- 
functions is complete in C e . 

4.2. Case (a) . In this Section we justify the leading terms of asymptotic expansions 
constructed above in case Ao ^ ob and thus v° ^ 0, see Section Let Ao be a 
solution to equation (|56"|) . All the functions {n, Nj, M, A4, V, 1Z, wq, ibi, v± and W2, 
V2) are as defined in Section FSTTI We also fix Ai according to l|57p. The approximate 
eigenvalue A e and eigenfunction W e are given by (fS"5)) and (fS"9")) respectively. 



, it does not satisfy the 



Notice that although W e £ if 1 (ft) since W e 

zero Dirichlet boundary conditions on dCl. To fix this we introduce the following 
boundary-layer corrector to our approximation. 



There exists a corrector V e solving the problem 
-div(a e W e ) + p e V e = in fl, 



Lemma 4.1 

(70) 

(71) V e \Ba = -W e \ga, V e 

such that U s = W e + V s G H^{fl) 



1 

and 



cm 

dn 



dV E 



dn 



max |y e I < Ce. 
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Proof. Clearly such solution of (f70|) . (|7T|) does exist. On each of the subsets Of and 
Qq the coefficients of (|70]l are smooth. Then the function 14 can reach its positive 
maximum or negative minimum only at the boundaries T e or dQ. Let us prove 
that this cannot be T £ . Suppose to the contrary the existence of x* <G T £ such that 
maxlVU = |V^(x*)|. The strong maximum principle yields that there is no more 

Q 

point inside {l\ or Q,q where the maximum is reached. Without loss of generality 
we assume V e (x*) > V e (x) for any x £ tt\T £ and V £ (x*) > (otherwise the point 
would be a positive maximum for — V £ and we would then consider — V £ ). Then by 
the virtue of Hopf 's Lemma [8, p. 330] applied in the relevant component of Qq we 
have 

f | <*-> >o - 

From transmission conditions (|7ip we have that the normal derivative on the f2f 
side of domain is also positive. Therefore the value of V £ increases from the point 
x* inside f2f in the n-direction and hence x* is not a point of maximum of V e in 
f2f. The contradiction proves that \V £ \ reaches it's maximum at dil. Then, from 
boundary conditions (f7Tj) , 



max I V E I = max I V £ I = max I W £ I < e max 
n an an on 



(x,J) 



e max 
on 



ii' i ( x, + eu> 2 ^ 



(-7) 



< Ce. 



Obviously U £ = W £ +V £ satisfies zero boundary condition on dfl and thus belongs 
to H^(Q). ' □ 

Lemma 4.2. The constructed corrector V £ satisfies the estimate \\V £ \\c^ < Ce 3 / 4 . 

Proof. Let X G C°°(R) and x{t) = 0, t < 1 and x(*) = 1, t > 2. Let us define a 
family of cut-off functions: 

Xs(x) = x (e- 1/2 dist(x,dQ)j , a; € ft. 

Then Xe '■ ft - * K satisfies the properties 

• Xe(x) = if dist (x, 9ft) < e 1 / 2 , 

• Xe{x) = 1 if dist (x, 90) > 2s 1 / 2 , 

• |Vx £ | < Ce" 1 / 2 and |suppVx e | < Ce 1 / 2 , 

where "supp" denotes a function support, and |supp-| is the measure of the corre- 
sponding support. Multiplying (|70[) by y 2 V^ and integrating by parts, we obtain 

(72) f a £ VV £ -V{ X 2 £ V £ )dx+ f p £X 2 V 2 dx = Q. 
Jn Jn 

Then using the identity 

w E • v(xlv £ ) = \V( X£ V £ )\ 2 - V 2 \V X e\ 2 , 

we get from (f72"j) 

(73) / a £ \V{ X eV £ )\ 2 dx+ f p £X 2 £ V 2 dx= f a £ V 2 \V X e\ 2 dx, 
Jn Jn Jn 

implying 

(74) f p £X 2 £ V 2 dx< f a £ V 2 \V X e\ 2 dx. 
Jn Jn 
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Lemma [4.11 provides the estimate V 2 < Ce 2 . Moreover, |suppVx e < Ce 1 / 2 and 
VXe| 2 < Ce~ 4 . Therefore estimate flU} yields 



(75) \\ X eV E \\h 



PeX. 



2 V 2 dx<Ce 3 ' 2 . 



Similarly we estimate 

(76) \\Q--Xs)Ve\\%= I p £ (l~Xe) 2 V 2 dx<Ce 3 / 2 , 

since |supp (1 - Xe)\ < Ce 1 ' 2 and \p e \ < Ce- 1 . Combining (JTSJ) and we obtain 

\\V e \\c = - Xe)V £ + XeVsWc* < CS^ 4 . 

Lemma 4.3. If p <E Hq(^) then 



□ 



(77) 



Jj^\ 2 dx) 1/2 < CMn 



Proof. Extend function ip by zero to whole of R" . Then (|77p follows upon rescaling 
y = x/s from the standard trace estimates applied to each connected component of 
Qo (which are shifts of Q ). □ 

Lemma 4.4. The corrected approximation U e satisfies the estimate 

\\B £ U £ - (A £ + l)- 1 ^!!^ < \\B £ U £ - (A £ + l^UeWn. < Ce 3 / 4 . 



Proof. For an arbitrary tp £ Ti. e consider 

\{BeU E - (Ae + l)" 1 ^,^! = lAe + lp'Kf/e - (A, + l)B e U e , ip) H .\ < 

<C\(U e ,cp)w-(Ae + l)(B e U e , l p)ne\=C\(U e , l p) n e-(Ae + l)(U e , <p)c* \ 



= C 
(78) < G 



a £ S7U £ V<p dx — A e / p e U e tpdx 
n Jn 



(div(a e V(7 e ) + A £ p e U e ) fdx 



+ C 



dUe 

dn 



dUe 

dn 



\ip\ dx. 



Denote the right-hand side of J7HJ) by F £ (U £ , tp). Substituting U £ = W £ + V £ and 
taking into account (|7D|) and (|7T|l. 



F £ (t/ £ ,^) <F £ (W/ £ ,<^)+C 



(A £ + 1) / Pe V e <pdx 



<F e (W c ,<p) + C\\V e \\c4<p\\ c . 



By Lemma 14.21 and obvious inequality ||<^||/;s < 

(79) F e (U e , ip) < F £ (W £ ,p) + Ce 3 ' 4 \\p\\ m - 

According to Lemma l3~2l (i) and (ii), 





8W e 




a 9We 




L 


E dn 





6 dn 


l 



F £ (W £ ,p) < Ce\\p\\ L 2 {n) +C 
Due to Lemmas 13.21 (Hi) and 14.31 the latter yields 
(80) F £ (W £ ,p) < Ce\\p\\ C t+Ce 3 / 2 (J^ \p\ 2 dx^' 2 < Ce\\p\\ H 

Using (|5D ]) . (fF9] l in ([75] ) yields 

\(B £ U £ - (A £ + l)- 1 ^, <p) w \ < Ce 3 / 4 \\p\\ w 
for all ip e H e . Hence, \\B £ U £ - (A £ + l)" 1 ^^ < Ce 3 / 4 . 



ip\ dx. 



□ 
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Lemma 4.5. \\U e \\c > C*^ 1 / 2 . 

Proof. By the triangle inequality we have 

(81) \\U e \\ w > \\U e \\ £ e > \\W e U.-\\V s \\ f 

We consider 

\\W e \\h=e 



1 J jwo (x, ^ + ewi (x, ^ + e 2 W2 (x, j dx+ 
+ J |v°(a:) + evi (x, — ^ + e 2 V2 (x, — ^ j da; = 
(82) = e _1 / Uo fa;, -)| 2 dx + 0(1) = e' 1 [ \v a (x)r) (-) f dx + O(l), e -> 0. 

./oe I V £/ I Joe I \ £ / I 



Extending : Qo — *• K by zero onto entire periodicity cell Q, by the mean value 
property we obtain 

(83) / v°(x)r] (-) 2 dx -> C* := {r] 2 ) [ \v°(x)\ 2 dx as e -> 0, 
where (f? 2 ) is the mean value of function rj 2 over Q, namely 

(V 2 ) = / V 2 {y)dy= / r l 2 {y)dy > 0. 

Since also v°(x) ^ 0, C* is positive. Therefore l|83")) and yield 

(84) ||W e |U. = CyV 1 / 2 + o^" 1 / 2 ), e^O. 

Due to LemmaiJand (21, it follows from flHTJ) that \\U e \\ C ' > C\ j2 s- 1 ! 2 +o(s- 1 l 2 ) 
as e — » 0. □ 

Theorem 4.6. Let A &e a solution to (|36p such that Xq ^ X® and Ai is defined 
according to (|57[) . TTien 

1. For sufficiently small e > £/iere exists an eigenvalue A e o/ @ smc/i that 

(85) |A £ -A -eAi| <Ci£ 5/4 , 

wf/j constant C\ independent of e. 

2. Let W e be defined by (|59[) and W e = ||W e ||7 e W e . XTien t/iere exist constants 
Cj(e) such that 

(86) W^-^C^uf <C 2 £ 5 / 4 , 

jeJ e £E 

where J e = {j : \Xj — Aq — eAj| < Ce 5 ^ 4 }, and Xj, Uj are eigenvalues and (C £ - 
normalized) eigenf unctions of and the constants C and Ci are independent of 
e. 

Proof. Application of classical lemma on "approximate eigenvalues" , e.g. [T8]. with 
U e = \\U £ \\^}U E as a test function and A e = A +eXi as an approximate eigenvalue, 
ensures, via Lemmas 14.41 & 14.51 the existence of an eigenvalue fx £ of operator B e 
such that 

(87) |(A e + l)- 1 -^1 <Ce b '\ 

and delivers the estimate analogous to (|86[) with u £ being eigenfunctions of B e and 

W e replaced by U £ . It suffices to notice that the eigenfunctions of the problem ([2J 
and of operator B £ coincide, their eigenvalues are related via = X £ + 1 and that 



16 



N.O. BABYCH, I.V. KAMOTSKI AND V.P. SMYSHLYAEV 



L e norm of the difference between U e and W e can be estimated via the right hand 
side of (JHU) (see Lemmas l4~2l and [475]) . □ 

Remark 2. Notice that (|86[) implies weaker but more transparent interpretations 
on the approximate eigenfunctions. For example, introducing 



u(x,y) 



v°(x), yeQ u 
w (x,y), y € Q , 



we claim that 
(89) 



(x,-) - Y d&X <Ce 3 / 4 



with appropriate dj(e). Note that \\u(-, -) ||i 2 (ri) > Co > 0. Then (|89[) follows 
from (|86[) by splitting its left hand side into the parts corresponding to £l\ and 
Qq, removing the weight, retaining only the main- order terms and then adding the 
inequalities up. 

We also remark that, in principle, the result on the convergence of eigenfunc- 
tions could be further sharpened, e.g. using the technique of two-scale convergence, 
cf. Section [5] below and [5]. 

4.3. Case (b). In this section we assume that Ao = Aj 3 for some j, its multiplicity 
is equal to 1 and the corresponding eigenfunction (j) has zero mean, i.e. (</>) = 0, see 
Section (3731 



Theorem 4.7. Let c € C 3 (f2), c = on dil, Ao be not a solution to (|36|) and 

Ai be defined according to (|B.17[) . Then there exist Eq > and constants C,C\ 
independent of £ (but dependent on c) such that for any < £ < £o, 
1. There exists an eigenvalue X e of @ such that 



(90) 



|A E -A -eAi| <Ce 5/4 



\W, 



Cj(e) such that 
(91) 



2. Let W* be defined by (JHSJ) and W' e 

^c 3 {e)u] e <Ci£ 5/4 , 



W* . Then there exist constants 



w; 



where J £ — {j : |A^ — Ao — eAx| < Ce 5 ^ 4 }, and A^,u|(x) are eigenvalues and (L e 
normalized) eigenfunctions of @. 



Proof. Proof of this theorem literally follows the proof of Theorem l4.6l with reference 
to Lemma \'S. 31 □ 



A direct analogue of Remark [2] also holds. 



5. On the eigenfunction convergence 

In this section we give a brief sketch of further refinement of the presented results 
using the technique of two-scale convergence, [TH [TJ Q15] . 

First, the inclusions intersecting or touching the boundary are "excluded", e.g. 
by re-defining a e and p e there as in the matrix phase (a e {x) — p £ {x) = 1). Denoting 



HOMOGENIZATION WITH DOUBLY HIGH CONTRASTS 



17 



now via e — > an appropriate subsequence in e, without relabelling, let u £ and A e 
be eigenfunctions and eigenvalues of the original problem, with normalization 



(92) / Vui + s z / Vui = 1. 

Jn\ Jo. e 

The boundedness of u £ in L 2 (fl) is then implied by (|92|) e.g. via the uniform 
positivity of the double-porosity operator whose form is given by the left hand 
side of (|92|) . [I9l Thm 8.1]. This implies that, up to a subsequence, w e — 1 u(x,y) 
and eVite — 1 V y u(x,y), where u € L 2 (Q, Hp er ) and — denotes weak two-scale 
convergence. Additionally, since (f9"2")l implies e||Vu e ||£2(n=) —> 0, [THJ, Thm 4.1] 
assures that the two-scale limit is independent of y in the matrix, i.e. is exactly in 
the form (88]). Further, by [H Thm 4.2], v a £ fl^(fi) and 

(93) 9\Vu e ^0 1 (y)(Vv°(x)+p(x,y)), 

where p £ L 2 (n, V pot ) with Q\ and 6\{y) denoting the characteristic functions of Of 
and Qi, respectively, and V pot denoting the space of potential vector fields on Qi, 
i.e. with respect to the Lebesgue measure supported on Qi, cf. [T^l §3.2]. 

Let X £ — > A and (X 6 — X )/e — > A x . Selecting then in ([2]) appropriate oscillating 
test functions <f> — <j) e one can pass to the limit recovering the weak forms of the 
equations derived in Section [3l For example, selecting 4> e {x) — eip(x)b(x/e), iji £ 
C °°(O), b(y) £ C™ r (Q) yields 

(Vv°(x)+p(x,y))-Vyb(y)'ilj(x)dydx+ / / \7 y wo(x,y)-V y b(y)ip(x)dydx = 



(94) = A / / w (x,y)b(y)tjj(x)dydx. 

Jn JQo 

This can be seen to be a weak form of (|23|) and (|2ip . Selecting further <fi £ (x) = ip(x) 
can be seen, after some careful technical analysis, to recover (15T1) . ([52)1 and (|54|) . 

The above implies that as long as (w ) 2 + ^ 0, Ao, Ai, v° and i^o can only be 
those constructed in Section [3] This does not however rule out the possibility that 
v° and wo are both trivial (equivalently, the two-scale limit u(x,y) is identically 
zero) . Therefore additional two-scale compactness type arguments are required, cf. 
[19l Lemma 8.2]. In fact, following literally the argument of Zhikov one observes 

that the two-scale compactness of the eigenfunctions does hold, i.e. u E — > u(x,y), 

2 

where — > denotes strong two-scale convergence, in particular there is a convergence 
of norms: 

(95) ||u s - u(x,x/e)||i2(n) as e -» 0. 

However, this in turn does not rule out the possibility of ||tt e || — > with the nor- 



malization (|92p . which requires a separate analysis. 

We announce here a partial result with this effect, postponing detailed discussions 
for future. 

Theorem 5.1. Let Xq be not an eigenvalue of the Dirichlet problem in Qq, i.e. 
Ao + Aj 3 , j > 1, see (|25J>. Then 

(i) In the above setting, necessarily, /3(Aq) > |Qi|Ao, i.e. there are gaps de- 
veloped for small enough e in the spectrum, containing in the limit at least 
{A : /3(A) < |Qi|A}. 
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(ii) If /3(Ao) = |Qi|Ao, necessarily u(x,y) ^ 0. Consequently, Ai can only be 
one of those described by (I57j) . The eigenf unctions converge strongly, in 
particular (|95j) holds. For fixed Xq and Ai, for small enough e the multiplic- 
ity of the eigenvalues X e near A £ = Ao + eAi coincides with the multiplicity 
of v as an eigenvalue of (JHH) , ([52"f . 

We remark that the above statement does not provide a full analogue of Hausdorff 
convergence of the spectra as in the double porosity case [TH Thm 8.1]. It does 
ensure however the existence of the gaps (on Figure [21 (Xf , Mj+i); 3 > 1) an d of 
the spectrum accumulation near the left ends /ij, j > 1, of the "bands" [/Xj-,A^]. 
However it does not clarify whether the "rests" of the bands, (fij, Aj 3 ] could be 
accumulation points. We conjecture that they could. For a chosen Ao = fij there 
exist infinitely many Ai = A^ according to (|57[) . (|56p . and X± — > +oo asn-> oo. 
On any band, for any small enough e there exists a finite but infinitely increasing 
number N(e) of eigenvalues according to ([85| . The issue is hence, in a sense, 
whether sX^ may become of order one for large n (n ~ N(e)). For A^ ~ e~ l , 
according to (|57| v ~ and hence, formally, the solutions v° of the homogenized 
equation ([51]) becomes oscillatory on the scale x/e 1 / 2 . One can attempt deriving 
asymptotic expansions similarly to those in Section 3, involving this new scale. A 
preliminary analysis has shown that those have formal solutions near every point 
inside the band. More detailed analysis is beyond the scope of the present work. 

Appendix A. Derivation of the limit equation for vq. 

Since 



— dy - - [ —n dy - [ d 2 fi 
dn x J r dxj 3 J Ql dxjdyj 



501) transforms to 



Taking into account (|4ip and ([3"5]) this becomes 



d 2 v f ON, 



dxjdxk J Ql dy 



(A x v + X v )\Q 1 \ + ^- I ^dy 



V3 



(A.1) =!?(-/ d ^y + [ V n jd y) + [ d ^ d y. 



dx i V Jqi d Vj Jr > Jr 9 

Since r)(y) = 1 on T, 

(A. 2) / ijnj dy — nj dy = 0. 



jii 



u 



We introduce homogenized matrix A hom = (A^l m ) r f k=l by (J53]) . According to ([46]) 
we have 

(A - 3) / dn- dV= dx- ~d^ dy + V dn~ dy + XlV dn~ dy - 

Substituting (jA.2j) - (|A.3[) into (|A.1[) yields 

(A.4) - divA^V^ = v(Xi)v° +K. jir - in fi, 

oxj 
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with 

KAi)=A |Qi|-Ax l^/y- j v ^dy 

and 

(A.5) /Cj = / -T^dy- —± dy. 

jQt 9 Vi Jr on y 

Lemma A.l. ^(Ai) depends on Ai with a non-zero linear coefficient. 
Proof. We estimate the linear coefficient 



dn y 



Note that 



where (l49l) and (l26l) have been used. Thus 



u(X 1 )=CX 1 +X \Q 1 \- f ?^dy 
Jr on v 

with positive constant C (depending on the choice of Ao). 

Corollary 1. (i) If \o = then rj(y) = 1 and hence C = \Qo\ 
(ii) According to (|48[) we also have the representation 

(A.6) i/(Ai) =CA 1 +A f|Q 1 |+ / PdyY 

Lemma A. 2. All K.j defined by (|A.5p equal zero. 
Proof. First we prove an auxiliary identity, namely 

dMj ? 97? 
<9n,, ' J dn. 



(A.7) / [^_^],,, / = ( , 

Notice for this that the left-hand side is 



dry 2 

-^—dy, 



I (AM jV - MjArj) dy = - f 2 1% dy = - 
where equations (|47[) and (f2"rj)) have been used. Since n(y) — 1 on T, 

/ 7T-dy=[ tf n i dy= [ rij dy = 0, 
•/Qo °Vi Jt Jr 

which finishes the proof of (|A.7[) . 
Then consider 

ICj = - I Nn dy- I ——2-T)dy, 



dn y 



which, according to P2"| and (|A.7[) . yields 
(A.8) Kj= I M^dy- f M 3 p-dy 



dn y a Jy 3 dn y 



Since Af and Nj are both harmonic in Qq, 
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Using (gT]) and (J43J) we obtain 



(A.10) 



dn 



Mj -r — dy = Nj - — dy 



dn 



a/v 



dn,. 



Substitution of (|A.9|I and (|A.10p into (|A.8p proves the lemma. 



□ 



Finally we come to the formulation of homogenized problem for the function vq, 
which comes from (|A.4[) and boundary condition (J4j) , resulting in |5T|) -(f52l). 

Appendix B. Proof of Lemma lcT3l 
We look for V\ , w\ , v 2 , W2 in the form 
(B.l) v 1 (x,y) = c(x)V 1 (y), 



(B.2) 
(B.3) 
(B.4) 



Wl (x,y) = c{x)W 1 {y) + ^-Z{ k \y), 

v 2 (x,y)=c(x)V 2 (y) + ^^V k (y), 
ax k 

w 2 (x,y) = c(x)W 2 (y) + ^-Z^(y), 

ox k 



where c is an arbitrary smooth function in f2 and Vi, Wi, z\ k ^ , Vk, i = 1, 2 , k 
1, ..,n, are functions to be found. 

Applying differential operator A e given by (f60|) to (|63|l in Of we obtain 

div (a £ VW*) + A £ p e W e * = 



= ^ e^cAyVi + E ° ( cA y V 2 (y) + j^- { A y T k (y) + ^}) + 



dyk 



(B.5) + e 1 ( 2 



d 2 v 2 



+ A x vi + Ao«i ) + s (A a: ii2 + Aiui + A0U2) + e A1U2 



Applying next A £ to ([63]) in Og we obtain 
div(a e VW*)+A e peW* = 

= |e (c [(A v + \ )Wi + \x<f>] + ^[(A, + Ao)^ fc) + 2^]) 



/ . , . _ d wi . . 

+ e ( (A y + Xo)w 2 + 2- — h A x w + XiWi 

OXjOyj 



(B.6) 



+ e z 2 



»2 



dxjdyj 



+ A x wi + X1W2 + e A X W2 



Evaluating the jumps of conormal derivatives on r e , we obtain 



dn 



dW F * 



On 



+ E 1 C 



84, dVi 
dn v dn v J 

(fe) 



jam _ av 2 \ dc l dz\ 



(B.7) 



+ E 



\ dn y dn y J ctefc | 9n 

^ 9«72 dV2 



+ n k c, 



dn y 



n k Vi 



Uer" 



+ 



V<9% dn* dn.y Nff ' ~ dn x l«fr 



+ £ 



3 9l02 



On the other hand function W e * is required to be continuous, i.e we have 
(B.8) Wi = Vi on r, 
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(B.9) Z (k) =0 on r, 

(B.10) W 2 = V 2 on I\ 

(B.ll) Z { 2 h) = V k on r. 

Equating to zero the term of order e _1 in (|B.5j) and the term of order e° in (|B.7I) . 
we obtain problem for Vi: 

(B.12) A„Vi=0 in Q X) |^ = |^ on T. 

(Jlty CfTty 

A solution to this problem exists since (0) = and we can present it as: 
(B.13) Vx=Vx+A, 

where < Vi >= and A is a constant which will be determined later. 

Equating to zero the term of order e° in (|B.6|) , and using (|B.8[) , (|B.9|) we obtain 
problems for Z± 



(B.14) (A y + A )Z} fc) =-2^- in Q , Z\ k} = on r, 



which admits an explicit solution 

(B.15) Z (k \y)=-y k <j>{y), 

and for W\ 

(B.16) (A y + A )Wi = -Ai0 in Q , Wi = Vi on T. 

A solution to the latter exists if and only if 

(B.17) A 1= / V 1 p-dy = - ( \WV 1 \ 2 dy = - f \VVi\ 2 dy, 

Jr dn y J Ql J Ql 

and we can present it in the following way: 
(B.18) Wi=m+Arj, 

where Wi solves problem (|B.16[) with V\ replaced by V\ (a solution exists for the 
same reason), and 77 solves (|26p (a solution exists since (<p) — 0). Notice that 

Mv) = - [ |^-rfy^o, 

Jr on y 



otherwise Ao would be a solution of (|36| which contradicts to the assumptions of 
this section. 

Equating to zero the term of order e° in (|B.5|) and the term of order e 1 in (|B.7[> . 
we obtain problems for V2 and Vk- For V2 we have: 

(B.19) A y V 2 = in Q x , —i = — ^ on V. 

OTly Oily 

A solution to this problem exists if and only if 

(B.20) 0= f^l da= ffl ds + Af^ d „, 

Jr on y Jr on y Jr on y 

and consequently 

On 



(B.21) A=(Ao(r ? »- 1 / ^dy. 
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The problem for V k has the form: 

ail pyn BZ^' 

(B.22) A y P k (y) = -2—± in Q u J-* = — L_ - n k V x on T. 

Oyk OTly UTly 

Solvability condition for this problem has the form 

The left hand side of ljR23]) can be transformed as follows, 
(B.24) / 2^ dy = -2 [ n fe Vi dy. 



Iq 1 dy k 

On the other hand, for the right hand side of (|B.23|) . 

' ' dzf 1 \ f ( d<t> 

nfeVi dy = / — 2/fco n kVi dy 



dn y I Jr \ dn y 

(B.25) = jf (-Vk^- - n k V^j dy = -2^ n k V x dy. 

Here we used (|B . 15|) . (IB.22[) and the integration by parts. Comparing (|B.24p and 
(|B.25|) we see that solvability condition (|B.23[) is satisfied. Finally W2 and Z% are 
arbitrary smooth functions satisfying (|B.10[) and (|B. 1 1[) . □ 
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